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The concept of parton fragmentation in QCD hard scattering phe- 
nomenology as well as NLO pQCD analysis of fragmentation functions are 
outlined. Hadroproduction of pions of a few GeV p± is discussed through 
the example of recent measurements at yj 5rhic = 200 GeV. 

1. Introduction 

Hard pQCD reactions produce highly virtual partons that will radiate off 
their virtualities, thus evolving into the non-perturbative states that contain 
observable hadrons. A quantitative understanding of this process, known 
broadly as the fragmentation of partons, considerably widens the class of 
reactions that can be handled within a perturbative QCD approach. Exam- 
ples are too many to enumerate here (as the reader may convince himself 
by browsing through the proceedings), two recent cases will be discussed 
below in Sections 3.2 and 3.3. 

A graph such as in Fig. 1 serves as a popular illustration of QCD fac- 
torization [1]. It also displays the role of fragmentation functions D^(z) in 
hard scattering phenomenology when there are detected hadrons in the final 
state. The early Field & Feynman [2] cascade model for quark fragmenta- 
tion into mesons is depicted in Fig. 2 side by side with a later field-theoretic 
definition in terms of a bilocal operator [3]. At leading twist accuracy, these 
describe fragmentation (or decay) functions (now in a narrower sense of the 
word) that turn a single parton into a single-hadron inclusive state, with 
fractional momentum transfer z from the parton to the hadron. Analogous 
to parton densities, renormalization induces evolution of the (transposed) 
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Fig. 2. Field & Feynman model (left) of cascade fragmentation and bilocal operator 
(right) a la Collins & Soper. 



DGLAP type: 

j=q,q,g 

This write-up discusses a few theoretical aspects of fragmentation func- 
tions and fragmentation processes. Please confer also K. Hamacher's con- 
tribution to these proceedings [4] for experimental results. 



2. Global Analyses of Fragmentation Functions 

Fragmentation functions are determined in NLO QCD fits to data. They 
are statistically dominated by e + e~ annihilations with particularly precise 
e + e~ —> hX data from Z° decays at LEP and SLD. Even though the e + e~ 
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data can be reproduced very successfully over a wide range of energy, there 
are several limitations: 

1. The leading order partonic process e + e~ —> qq produces quarks only; 
the subleading terms are too weak to constrain gluon fragmentation 
sufficiently. Gluon-quark mixing in the evolution Eq. (1) is not enough 
of a constraint, either. 

2. Data are less accurate at large-z. This fact translates into less accurate 
functions A(^)li arge -^ 

3. The individual light flavour channels (corresponding to the "Ranks" 
in Fig. 2) cannot be disentangled from the flavour inclusive data. 

These items limit, at present, the accuracy of the functions D{(z). They re- 
quire that fragmentation analyses become truly global the way that analyses 
of parton densities already are, i.e. including data from different processes 
that constrain different kinematic regions and flavour combinations. 

Recent global analyses of fragmentation functions, determined at Next- 
to-Leading-Order accuracy in QCD, were performed in [5]. The individual 
fits are conceptually very similar, such that a comparison between the fits 
gives a good first estimate of the uncertainties in the functions Di(z). As 
explained above, the uncertainties of these functions are data driven and can 
be large. It can indeed be shown [6] that the electroweak couplings of up 
and down type quarks at the Z° pole suggest that only flavour-insensitive 
combinations such as the singlet J2 q D q +q are well determined. Flavour 
structure can be constraint through the analysis of semi-inclusive DIS [6]. 
The simple lowest order term 

c^sidis oc W) D ^ Z ) + Q(x)D q (z)] (2) 
q 

illustrates that the information is local in z (no convolution at LO) and that 
the by now very good knowledge of the parton distributions q(x) provides 
reliable flavour weights. 

Having discussed item 3, we are still left with items 1, 2 from our shop- 
ping list above. They lead over to the analysis of hadroproduction reactions. 

3. Fragmentation in Hadroproduction Reactions 

3.1. Partonometry at VSrhic = 200 GeV 

RHIC (pp mode) operates at \J Srhic — 200 GeV, between lower energies 
(such as ISR or Fermilab), where the pQCD calculations are known to be 
problematic [7], and higher CMS energies (up to V^Tevatron), where pQCD 
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Fig. 3. Partonic decomposition of the initial and final state of pp collisions at 200 
GeV; central and forward rapidity. 
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Fig. 4. Scaling variables, averaged over the parton initial and final state of pp 
collisions at 200 GeV; central and forward rapidity. 



seems to work [5]. It was quite remarkable to see, how well pion production 
data in RHIC pp collisions [8] agree with the NLO QCD calculation [9] 
at central and forward rapidity and down to p± values as low as p± > 
1 GeV. While it remains a challenge to understand the y/~S dependence 
of hadroproduction spectra in terms of power corrections [10], the good 
agreement of the pion spectrum with NLO QCD at \J Srhic and the good 
precision of the RHIC data are a goldmine for fragmentation analyses. 

Fig. 3 decomposes RHIC collisions into 2-parton initial and 1-parton- 
inclusive final states. For Fig. 4 the corresponding scaling variables have 
been evaluated at their average values, manifesting the symmetric and asym- 
metric kinematics of central and forward rapidities, respectively. The parton 



kretzer printed on February 2, 2008 



5 



densities are probed in regions of x\^ where they are well known already 
whereas new constraints are imposed on the large- z fragmentation func- 
tions, and in particular on the gluon function. While the constraint on the 
gluon fragmentation function is nonlocal (convoluted) in hadroproduction 
reactions, a QCD analysis of this reaction is more straightforward than an 
analysis of anti-tagged gluon jets in e + e~ — ► bbg configurations [4]. The 
latter are nontrivial due to the fact that jets and fragmentation functions 
are not the same in QCD. 

I will next leave the fragmentation functions and move on to discussing 
theoretical aspects of two recent pion production measurements at RHIC. 



3.2. The double-spin asymmetry A^ L 

The gluon helicity contribution Ag to the proton spin is probed in colli- 
sions of longitudinally polarized protons under the condition that the domi- 
nant dynamics are perturbative and of leading twist origin. These conditions 
may apply to the measurement of beam helicity asymmetries 

dAa da ++ — da^ 
= ~do~ = da+++da+- (3) 

for high transverse momentum reactions where da = da / dp±\ l8iYge _ p± . For 
large p± pion production, a variant of Fig. 1 with polarized beams illus- 
trates that the (glue- glue contribution to the) reaction can schematically be 
written as 

A LL oc ( Ag a Ag h ) Aa ab ^ cX <g> . (4) 
The Mellin transform 

Ag a (N)Ag b (N)a ab ^ cX (N)D«(N) (5) 

[with Ag a (N) = x N ~ 1 g a (x) etc.] turns the nonlocal (in parton momen- 
tum, i.e. x a 7^ Xb) convolutions 

Ag a ® g a • • • = / dx a dx b Ag a (x a )Ag b (x b ) ... (6) 

a 

into products that are local in the moments N. For small to medium p± < 4 
GeV one can then show [11] by variation SA^/SAg = that A^ L has a 
minimum and is basically positive definite at leading twist, i.e. 

^lI p± <4GcV><?(-10- 3 ) (7) 

is a necessary condition for leading twist dynamics to be dominant in the 
production of small to medium p± pions. To achieve a slightly negative AJ L 
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of order O ( — 10 -3 ) the Mellin inversion of Ag a ^(N) back to Ag a ^(x a ^) has 
to generate a crossing of Ag{x) through zero and explore the fact that the 
convolution in Eq. (6) can then pick x a ^ x^ such that Ag a {x a )Agi ) {xi ) ) < 0. 
For any foreseeable experimental accuracy O (— 10 -3 ) ~ which is why 
A£ L is "basically" positive definite. At present, data [12] are statistically 
not conclusive and the positivity of A^ L remains yet to be confirmed, or 
disproved. 

3.3. Production of p(— > 2tv) mesons 

Recent data [13] indicate that the invariant mass spectrum of pion pairs 
in pp collisions peaks at a value around, but significantly below, the p mass 
as measured in e + e~ exclusive production. Now, resonance mass shifts 
would seem an unlikely manifestation of a perturbative effect if it were not 
for the fact that the effect extends to large values of transverse momen- 
tum p± ^ 3 GeV where the p can be expected to be a parton fragment. 
The strong decay p — > 2tt then contributes a p-wave resonance to the D? n 
fragmentation function [14]. 1 

In Fig. 1, the fragmenting parton momentum k runs over the following 
virtualities and transverse momenta: 

rs r(l-z)/z(k 2 -M 2 /z) 

/ dk 2 / dk\ . (8) 

Jm 2 /z Jo 

Strictly, Eq. (8) differs from the operator product depicted in Fig. 2 which 
is defined to be UV-divergent [after integration over J d 4 kS(k + — P + / z)]. 
The argument I am about to make, though, is not related to the subtle [15] 
problems that arise from cutting off momentum integrals. 

The partonic center of mass energy s can, for an order of magnitude 
estimate, be inferred from Fig. 4 through (s) = {x\X2) Srhic- Under the 
assumption that the partonic spectrum behaves like A: -2 , a Breit-Wigner 
peak in the 2tt mass spectrum - positioned at the p mass for p± —> oo - is 
then shifted to a lower value by a few 10 MeV around a finite p± = 3 GeV. 
This suggests a partonic effect that is dual to the pion phase space distortion 
at low p± —> 0. A more detailed evaluation will be given in [16]. 

4. Conclusions 

Parton fragmentation processes explore conceptually and phenomeno- 
logically rich aspects of hadronization, as I hope to have demonstrated 

1 For the semi- quantitative argument we are making here, we will ignore the mixing 
of m-hadron fragmentation with n-hadron (n < m) fragmentation, as well as the 
non-resonant s wave channel. 
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through the examples of a positivity bound on the spin asymmetry AJ L 
and a perturbative ansatz to explain resonance mass shifts at high p±. I 
have also outlined the next steps in global analysis of fragmentation func- 
tions; a corresponding update will be made available very soon. 

5. Acknowledgments 

Many thanks to M. Grosse-Perdekamp for the invitation and to all orga- 
nizers of ISMD04, in particular B. Gary, for financial support and a superbly 
organized workshop. This write-up benefited from discussions and corre- 
spondences with A. Bacchetta, K. Barish, S. Brodsky, G. Bunce, P. Fachini, 
K. Hamacher, E. Laenen, A. Mueller, M. Radici, E. Sichtermann, R. Seto, 
and W. Vogelsang. Work supported by RIKEN-BNL and the U.S. Depart- 
ment of Energy (contract DE-AC02-98CH10886). 

REFERENCES 

[1] J.C. Collins, D. Soper and G. Sterman, in A.H. Mueller, ed., Perturbative 
Quantum Chromodynamics (World Scientific 1989), hep-ph/0409313. 

[2] R.D. Field and R.P. Feynman, Nucl. Phys. B136, 1 (1978). 

[3] J.C. Collins and D.E. Soper, Nucl. Phys. B193, 381 (1981), B213, 545 (1983) 
(E); B194 445 (1982). 

[4] K. Hamacher, these proceedings. 

[5] S. Kretzer, Phys. Rev. D62 054001 (2000); B. Kniehl, G. Kramer and 
B. Potter, Nucl. Phys. B582 514 (2000); Nucl. Phys. B597 337 (2001); 
L. Bourhis, M. Fontannaz, J.P Guillet, M. Werlen, Eur. Phys. J. C19 89 
(2001). 

[6] S. Kretzer, E. Leader, E. Christova, Eur. Phys. J. C22, 269 (2001). 
[7] C. Bourrely and J. Soffer, hep-ph/0311110. 
[8] E. Sichtermann, these proceedings. 

[9] B. Jager, A. Schafer, M. Stratmann and W. Vogelsang, Phys. Rev. D67, 
054005 (2003); and references given there. 

[10] G. Sterman, these proceedings; G. Sterman and W. Vogelsang, 
hep-ph/0409234; D.W. Sivers, S.J. Brodsky and R. Blankenbecler, 
Phys. Rept. 23, 1 (1976). 

[11] B. Jager, M. Stratmann, S. Kretzer, W. Vogelsang, Phys. Rev. Lett. 92, 
121803 (2004). 

[12] PHENIX Collab. (S.S. Adler et al), hep-ex/0404027. 

[13] STAR Collab. (J. Adams et al), Phys. Rev. Lett. 92, 092301 (2004); and 
references to theoretical work given there. 

[14] A. Bacchetta and M. Radici, hep-ph/0409174; and references given there. 



8 



KRETZER PRINTED ON FEBRUARY 2, 2008 



[15] J.C. Collins, Acta Phys. Polon. B34, 3103 (2003). 

[16] A. Bacchetta, A. Dumitru, S. Kretzer, and M. Radici, in preparation. 



